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Abstract 



The positions of columnar pins and magnetic flux lines determined from a 
decoration experiment on BSCCO were used to calculate the single-particle 
density of states at low temperatures in the Bose glass phase. A wide Coulomb 
gap is found, with gap exponent s ~ 1.2, as a result of the long-range in- 
teraction between the vortices. As a consequence, the variable-range hop- 
ping transport of flux lines is considerably reduced with respect to the non- 
interacting case, the effective Mott exponent being enhanced from po = 1/3 

to Pcflf ^ 0.5 for this specific experiment. 
PACS numbers: 74.60.Ge, 05.60.-hw 
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The remarkably rich phase diagram of magnetic flux lines in high-temperature super- 
conductors, especially when subject to point and/or extended pinning centers, has attracted 
considerable experimental and theoretical interest [jl|. Understanding the interaction of 
vortices with defects is especially important since flux lines must be pinned to mimimize 
dissipative losses from flux creep. A promising pinning strategy involves the creation of lin- 
ear damage tracks in materials by heavy ion irradiation. These columnar defects have been 
found to increase the critical current, as well as to shift the irreversibility line significantly 
upwards 0. 

A theory of flux pinning by correlated disorder has been developed to explain these re- 
sults, exploiting a formal mapping of the statistical mechanics of directed lines onto the 
quantum mechanics of two-dimensional bosons 0. In this study, the intervortex repulsion, 
whose range is determined by the London penetration depth A, was only treated using ap- 
proximate, order of magnitude estimates. However, if A ^ oq, where = (4/3)^/'^(0o/-B)^/^ 
is the average distance between vortices (0o = hcjle is the elementary flux quantum), the 
interactions become effectively long-range, and may lead to important correlation effects. 

Indeed, the analogy of flux lines at low temperatures pinned to columnar defects (Bose 
glass phase) with two-dimensional localized carriers in doped semiconductors (Coulomb 
glass) 0, suggests that a "Coulomb" gap should emerge in the single-particle density of 
states (distribution of vortex pinning energies). Because such a gap will affect significantly 
the current- volt age characteristics in a variable range hopping approach , it is important 
to estimate the size of the Coulomb gap in the Bose glass phase and to understand the 
correlation effects induced by the intervortex repulsion. The recent successes, moreover, 
in simultaneously measuring both the columnar defect and flux line positions [^,0 allow 
for detailed comparison of the spatial correlations found in experiment and theory. Such a 
comparison is not feasible in the semiconductor case. 

The density of states may be obtained by using a variant of the Monte Carlo algorithm 
described by Shklovskii and Efros [^0. Using the experimentally determined columnar 
defect and flux line positions shown in Fig. 1, we can predict the density of states and the 
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transport characteristics for this specific sample in the variable range hopping regime, at 
temperatures slightly below the depinning transition. We find that the ensuing Coulomb 
gap is remarkably large, even in the case A ~ ao, i.e. when the interactions are relatively 
short-range. Vortex interactions raise the effective Mott exponent from the non-interacting 
result Pq = 1/3 to a value p^s ~ 0.5 for this specific sample. These results imply that 
correlation effects strongly enhance the pinning of flux lines to columnar defects. 

We consider the following model free energy for flux lines, deflned by their trajectories 
ri{z) as they traverse the sample of thickness L, with the magnetic fleld B aligned along the 
z axis (perpendicular to the CuO planes), 



Here rij{z) = \ri{z) — rj{z)\, and V{r) = 2eoKQ{r / X) is the repulsive interaction potential 
between the lines; the modifled Bessel function Kq(x) oc — Inx for x — 0, and oc x^^^^e"^ 
for X — >■ oo. Thus the (in-plane) London penetration depth A deflnes the interaction range. 
The energy scale is set by eo = (0o/47rA)^, and ei is the tilt modulus. Finally, the pinning 
energy is a sum of Nd z-independent potential wells Vn with average spacing d centered 
on the {Rfc}, whose diameters are typically Cq ~ lOOA, with a variation induced by the 
dispersion of the ion beam of Sck/cQ ~ 15%. This induces some probability distribution 
P of the pinning energies Uk, which may be determined from the (interpolation) formula 
Uk ~ (eo/2) ln[l + (cfc/\/2^)^], where ^ is the coherence length [|]. E.g., for A ^ 4200A one 



formally mapped onto a two-dimensional zero-temperature quantum mechanical problem 
by using a transfer matrix approach. In this boson analogy, the real temperature T assumes 
the role of h in the quantum problem, and the boson electric fleld and current density map 
on the superconducting current J and the true electric fleld £, respectively (see Table I in 




(1) 




As is explained in detail in Ref. 0], the statistical mechanics of the model (^ can be 



3 



Ref. 0). Thus the roles of conductivity and resistivity become interchanged. 

We are interested in the low-temperature properties of flux lines pinned to columnar 
defects, with filling fraction / = N/No = (d/ao)'^ < 1, in the Bose glass phase, where 
all the vortices are assumed (and found) to be bound to the pinning centers. For T less 
than a characteristic fiuctuation temperature Ti, one arrives at the classical limit of the 
corresponding boson problem (/i — > 0), and as the fiux lines are now well separated, the Bose 
statistics become irrelevant. Furthermore, in this limit thermal wandering is suppressed, and 
the fiux lines will be essentially straight; hence the tilt energy in Eq. (JlD can be neglected. 
In the boson representation we eventually have to deal with a time-independent problem 
defined by the two-dimensional effective Hamiltonian 

]^ Nd Nd 

H = -J2 nirijVinj) + niU , (2) 

where i,j = 1,...,Nd denote the defect sites, randomly distributed on the xy plane. 
rii = 0, 1 is the corresponding site occupation number, and (originating in the varying 
pin diameters) the tj are random site energies, whose distribution may be chosen to be cen- 
tered at t = 0, with width w [for simplicity, we assume a flat distribution P{ti) = l/2w for 
\ti\ < w, and P{ti) = else]. 

The Hamiltonian (H) is precisely of the form studied in the context of charge carriers lo- 
calized at random impurities in doped semiconductors (Coulomb glass problem) @J^,|]. We 
remark that it is equivalent to a two-dimensional random-site, random-fleld antiferromag- 
netic Ising model with long-range exchange interactions, and has, at least to our knowledge, 
eluded successful analytical approaches going beyond simplifying mean-fleld type consid- 
erations [^], and phenomenological scaling arguments ||^. Therefore one has to resort to 
numerical studies using suitable Monte Carlo algorithms, as described in Refs. 0,0. 

Basically, starting from a given distribution of defect sites, N < of which are oc- 
cupied, single-particle energies are calculated according to e, = dH/dui = Ylj^i njV{rij)+ti. 
The initial conflguration is then relaxed by moving single "particles" to empty places, until 
precisely the N lowest energy levels are occupied (eJnax = niax„.=iei < e^^j^ = min„_oej)- The 
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ensuing intermediate state is then probed against all possible single-particle hops from filled 
to empty defect sites, with associated energy change Ai_j = ej—ei — V{rij). If any Aj^^ < 0, 
the move from site i to j is performed, and thus the total energy is decreased. (Note that this 
procedure constitutes a zero-temperature algorithm.) Afterwards the "equilibration" step 
has to be repeated, for all the of course rearrange upon changes in the occupation numbers 
rii. Finally, a chemical potential /i is calculated (approximately) by /i = {e^^^ — emax)/^- 
[Formally, this corresponds to adding a term —fij^i^i to Eq. (H); this chemical potential is 
related to the external magnetic field H hj fi = eoln(A/^) — if0o/4vr — (Uk), where (Uk) 
includes a small thermal renormalization of pinning energies.] 

By repeating this procedure for different initial configurations, one may then obtain the 
density of states g{e) from the site energy statistics. We shall use a normalization of the 
density of states such that / g{t)dt = l/dP. To find the correct ground state for each con- 
figuration, one would in principle have to test the configuration against any simultaneous 
ra-particle hops, n = 2, . . . , oo. However, previous investigations have shown that termi- 
nating at n = 1 already yields at least qualitatively reliable estimates for the energy level 
distributions We have performed many such simulations for a variety of filling frac- 

tions / and values of X/d. Details of these investigations and their results shall be reported 
elsewhere [|lOl. In this Letter, we shall instead use these techniques to determine the density 
of states directly from experimental data, obtained as follows. 

The positions of columnar defects and fiux lines in a Bi2Sr2CaCu208 (BSCCO) crystal 
were determined simultaneously by a combined chemical etching/magnetic decoration ap- 
proach (for details, see Ref. 0). A digitized SEM image of the positions of = 162 fiux 
lines and Nd = 686 columnar defects (hence / ~ 0.24) for a BSCCO sample irradiated 
with corresponding matching field = 118G {d ^ 4400A) and decoration field B = 27G 
(ao ~ 9400 A) is shown in Fig. 1. The critical and irreversibility temperature under these 
conditions are Tc = 87K and = 81K, respectively. Assuming the fiux line distribution is 
frozen in at Tin-, we estimate that the effective London penetration depth is A(Tii.r) ~ 4200 A 



TT[, and X/d ^ 0.96. 



In Fig. 2 the structure factor S{q) for the flux hues is depicted, as obtained from S{q) — 
averaging over directions in Fourier space (thick hne) . Clearly the flux line 
distribution is highly correlated, with S{q) displaying a peak at go^o ~ 27r. We have used the 
experimental defect and flux line positions as an initial state for the Monte Carlo routine. 
In order to minimize boundary effects, we have kept the configuration fixed in a frame 
extending 10% (which amounts to ~ 3A) inwards from each of the rectangular boundaries. 
This leaves 464 sites and 106 displacable particles for the simulation. For vanishing site 
randomness (w = 0), it turns out that about 20% of the flux lines are moved in the course of 
the energy minimization process, which may in part be attributed to finite size effects, and 
also to the unaccounted variation in sample pinning energies. Using the more realistic value 
w — O.leo instead, we find typically 40% changes with respect to the original experimental 
distribution. But, as can be seen in Fig. 2, in both cases the highly correlated character of 
the flux line distribution is preserved, although the height of the peak at qo decreases upon 
increasing w. We have found, however, that substantially stronger disorder would destroy 
the peak in S{q). Similarly, for considerably lower values of X/d {X/d < 0.2), the spatial 
correlations also disappear, because the random site energies would then dominate over the 
interactions. The correlations we do find strongly support the assumption that the effective 
London penetration depth is the larger value Aeg = A(T = Tin). 

The density of states for the 464 "inner" sites, as obtained from averaging over 100 runs 
with different assignments of random site energies, drawn though from the same distribution 
P{ti) with width w = O.leo, and using Xes for the interaction range, is shown in Fig. 3. 
(Results for the density of states are remarkably insensitive to the precise value w.) The 
Coulomb gap, separating the occupied and empty energy levels, is remarkably wide, its 
width amounting to ~ O.Scq at half maximum [d'^g[e) ~ 0.8], i.e., almost a third of the total 
width of g{e) there. Near its minimum, this pseudogap may be described by the formula 

g{e) (X \e - , (3) 

with a gap exponent s fa 1.2 for the specific parameter values here. In the immediate vicinity 
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of the chemical potential, this power law is smeared out, and g{fi) is actually finite (but very 
small) due to the finite range of the interactions. We remark that in the limit A — > oo and 
for small filling, s ~ 3 may be reached . 

From this single-particle density of states, we may now infer the transport properties 
in the variable-range hopping regime by minimizing the free energy of two superkinks of 
size R and separation Z, 6Fsk = 2EKR/d + ZA{R) — /lRZ [||. Here, the first term 
consists of the energy of the double kink, with Ek = y/iiUod, and the third one derives 
from the Lorentz force = ipoJ/c induced by the current J. The second contribution 
stems from the fact that for a hop of distance i?, the available energy states lie in the 
interval [/i, A(i?)], where in D dimensions (here D = 2) A{R) is determined by the equa- 
tion g{e)de = R~^ . Minimizing first for J = gives the longitudinal extent of the 
kink to be Z* = —2Ek / d{d/S./ dR)R* . To first order in J, one subsequently arrives at 
J0o/c = A{R*)/R*, and thus 6F* = 2EKR*{J)/d is the result for the optimized free en- 
ergy. The latter finally enters the resistivity p as an energy barrier in an activation factor, 
S = PqJ exp{—6F* /ksT). Using g{e) of Fig. 3 then amounts to studying transport slightly 
below the depinning temperature, at T = 77K, say. Note that for this field range ii ~ eo and 
thermal renormalizations of pinning energies become relevant only for Ti ^ 0.9Tc (~ 78K 
here) [3] . Thus, all temperatures T < Ti may be considered as "low" . 

In the regime where Eq. (^) holds, the final result from these considerations for the highly 
nonlinear current- volt age characteristics may be cast in the form 

S ^ poJexp [-{2EK/kBT){Jo/jy] , (4) 

where p is an exponent generalizing Mott's law [pq = 1/{D -\- 1)], which is valid in the 
case of vanishing interactions. For long-range interactions producing a Coulomb gap of the 
form (Q), one finds p = {s + 1)/{D -\- s + 1). Fig. 4 shows a log-log plot of the function 
R*{J)/d vs. j = J(t>od/2eoc, derived from the density of states in Fig. 3, as compared to a 
similarly calculated curve with the vortex repulsion being switched off. While in the latter 
case the result is indeed a straight line with slope —1/3, interactions considerably enhance 



the pinning by raising the effective Mott exponent to pes ~ 0.5 for not too low values of J. 
For J the cutoff of the interaction at A reduces p^s somewhat. 

Similar to the gap index s, from which it is derived, p really should not be understood as 
a universal number, but rather as some effective exponent pes conveniently describing the IV 
characteristics. Its value in general depends on both the filling / and the interaction range 
X/d; its maximum value pes ~ 0.68 is reached for A oo and small / [|lOl- These results 



clearly rule out the mean-field estimate, which would yield s = D/a — 1 and p = 1/(1 + 0") 
for a potential V{r) oc r^"' [a < D\ a. logarithmic interaction is recovered in the limit 
(J — >■ 0). Rather, they seem consistent with the analysis by Fisher, Tokuyasu, and Young, 
who relate p with the stiffness exponent ~ —0.5 and a fractal index 1 < ipq < 2 for 
the supposedly equivalent gauge glass model. Their scaling relation p = 1/(1 + \Q\/iPq), 
indeed yields p ^ 2/3 for ■^/'g ^ 1. It is, however, not yet settled if universality does apply in 
the purely long-range limit, and even the results by Mobius and Richter, who were able 
to simulate the 1/r Coulomb problem in two and three dimensions for very large systems, 
have apparently not reached the fully asymptotic regime. Also, our derivation of R*{J) from 
g{e) of course constitutes an approximation which neglects some subtle correlations, e.g. the 
spatial clustering of those sites which are energetically close to fj,, |]7|JTI|, and may be subject 
to corrections in the limit J — 0. 

In summary, we have demonstrated that the vortex-vortex repulsion can lead to re- 
markable correlations both in real space and in the single-particle density of states, even 
for A ~ Oq, which is easily accessible in experiment. An important consequence of these 
correlation effects is the drastic enhancement of flux line pinning to columnar defects in the 
Bose glass phase, whenever A > oq. 

We benefitted from discussions with A.L. Efros, D.S. Fisher, T. Hwa, P. Le Doussal, and 
V.M. Vinokur. This work was supported by the National Science Foundation, primarily by 
the MRSEC Program through Grant DMR-9400396, and through Grants DMR-9417047 and 
DMR-9306684. U.C.T. acknowledges support from the Deutsche Forschungsgemeinschaft 
(DFG) under Contract Ta. 177/1-1. 
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FIGURES 



FIG. 1. Positions of columnar pins (open circles) and flux lines (filled circles), as obtained 
from the experiment. 

FIG. 2. Flux line structure function S{q) as obtained from the experiment (thick line), 
and from the simulation with w — (dashed), and w — O.leo (thin line), averaged over 100 
different assignments of random site energies. 

FIG. 3. Normalized density of states G{E) = (fg{e) as function of the single-particle 
energies E — e/2eo, averaged over 100 runs {X/d — 0.96, w — O.leo). The mean chemical 
potential is ///2eo — 0.65. 

FIG. 4. Double-logarithmic plot of the exponential factor R*{J)/d for Mott variable- 
range hopping vs. j — J(j)od/2€oC, as obtained from the density of states in Fig.3 (filled 
circles) , compared to the result of the non-interacting case (open circles) . 
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Figure 1 
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